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Abstract 

O _ 

■ The charmonium and bottomonium mass spectra are calculated with the 
On . systematic account of all relativistic corrections of order v'^/c^ and the one- 
loop radiative corrections. Special attention is paid to the contribution of 

^ I the retardation effects to the spin-independent part of the quark-antiquark 

potential, and a general approach to accounting for retardation effects in 

■ the long-range (confining) part of the potential is presented. A good fit to 
available experimental data on the mass spectra is obtained. 
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I. INTRODUCTION 



p • The investigation of the meson properties in the framework of constituent quark models 

O , is an important problem of the elementary particle physics. At present a large amount of 
experimental data on the masses of ground and excited states of heavy and light mesons has 
been accumulated By comparing theoretical predictions with experimental data, one 
^ ■ can obtain a valuable information on the form of the quark-antiquark interaction potential. 
■ Such information is of great practical interest since at present it is not possible to obtain 
the qq potential in the whole range of distances from the basic principles of QCD. As it is 
well known, the growing of the strong coupling constant with distance makes perturbation 
theory inapplicable at large distances (in the infrared region). In this region it is necessary 
to account for nonperturbative effects connected with the complicated structure of the QCD 
vacuum. All this leads to a theoretical uncertainty in the qq potential at large and interme- 
diate distances. It is just in this region of large and intermediate distances that most of the 
basic meson characteristics are formed. This makes it possible to investigate the low-energy 
region of strong interaction by studying the mass spectra and decays of mesons. 

Some recent investigations have shown that there could be also a linear (in radius) 
correction to the perturbative Coulomb potential at small distances (in contradiction with 
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OPE predictions). The estimates of the slope yield that it could be of the same order of 
magnitude as the slope of the long-range confining linear potential. It means then that the 
widely used Cornell potential (the sum of the Coulomb and linear confining terms) is really 
a correct one in the static limit both at large and at small distances. 

The relativistic properties of the quark-antiquark interaction potential play an impor- 
tant role in analysing different static and dynamical characteristics of heavy mesons. The 
Lorentz-structure of the confining quark-antiquark interaction is of particular interest. In 
the literature there is no consent on this item. For a long time the scalar confining ker- 
nel has been considered to be the most appropriate one The main argument in favour 
of this choice is based on the nature of the heavy quark spin-orbit potential. The scalar 
potential gives a vanishing long-range magnetic contribution, which is in agreement with 
the flux tube picture of quark confinement of 0, and allows to get the fine structure for 
heavy quarkonia in accord with experimental data. However, the calculations of electroweak 
decay rates of heavy mesons with a scalar confining potential alone yield results which are 
in worse agreement with data than for a vector potential The radiative M 1-transitions 
in quarkonia such as e. g. J/ip ^ rjc'y are the most sensitive to the Lorentz-structure of 
the confining potential. The relativistic corrections for these decays arising from vector and 
scalar potentials have different signs J^,^. In particular, as it has been shown in ref. [Q, 
agreement with experiments for these decays can be achieved only for a mixture of vector 
and scalar potentials. In this context, it is worth remarking, that the recent study of the qq 
interaction in the Wilson loop approach indicates that it cannot be considered as simply 
a scalar. Moreover, the found structure of spin-independent relativistic corrections is not 
compatible with a scalar potential. A similar conclusion has been obtained in ref. |T^ on 
the basis of a Foldy-Wouthuysen reduction of the full Coulomb gauge Hamiltonian of QCD. 
There, the Lorentz-structure of the confinement has been found to be of vector nature. The 
scalar character of spin splittings in heavy quarkonia in this approach is dynamically gen- 
erated through the interaction with collective gluonic degrees of freedom. Thus we see that 
while the spin-dependent structure of (qq) interaction is well established now, the spin-inde- 
pendent part is still controversial in the literature. The uncertainty in the Lorentz-structure 
of the confining interaction complicates the account for retardation corrections since the 
relativistic reconstruction of the static confining potential is not unique. In our previous 
paper |1T| we gave some possible prescription of such reconstruction which, in particular. 



provides the fulfilment of the Barchielli, Brambilla, Prosperi (BBP) relations [|T2| following 
from the Lorentz invariance of the Wilson loop. Here we generalize this prescription and 
discuss its connection with the known quark potentials and the implications for the heavy 
quarkonium mass spectra. 

The other important point is the inclusion of radiative corrections in the perturbative 
part of the quark potential. There have been considerable progress in recent years and now 
the perturbative QCD corrections to the static potential are known up to two loops p!3|JT^ 
though for the velocity dependent and spin-dependent potentials only one-loop corrections 
are calculated [P"3|-P^. 

The paper is organized as follows. In Sec. II we describe our relativistic quark model. 
The approach to accounting for retardation effects in the qq potential in the general case is 
presented in Sec. HI. The resulting heavy quark potential containing both spin-independent 
and spin-dependent parts with the account of one-loop radiative corrections is given in 
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Sec. IV. We use this potential for the calculations of the heavy quarkonium mass spectra in 
Sec. V. Section VI contains our conclusions and discussion of the results. 



II. RELATIVISTIC QUARK MODEL 

In the quasipotential approach a meson is described by the wave function of the bound 
quark-antiquark state, which satisfies the quasipotential equation of the Schrodinger 
type 



\ J {2tt 

where the relativistic reduced mass is 



^m(p) = / j^,V{p,q;M)^M{q), (1) 



and Ea, Ef, are given by 



_ — m? + „ — ml + m? 

Here M = Ea + E^ is the meson mass, irta^b are the masses of light and heavy quarks, and p 
is their relative momentum. In the centre of mass system the relative momentum squared 
on mass shell reads 

The kernel V(p, q; M) in Eq. (|I|) is the quasipotential operator of the quark-antiquark in- 
teraction. It is constructed with the help of the off-mass-shell scattering amplitude, projected 
onto the positive energy states. Constructing the quasipotential of the quark-antiquark in- 
teraction we have assumed that the effective interaction is the sum of the usual one-gluon 
exchange term with the mixture of long-range vector and scalar linear confining potentials, 
where the vector confining potential contains the Pauli interaction. The quasipotential is 
then defined by 

+Vy{y)V^j:,,^ + l^5(k)|M„(g)^Zb(-g), (5) 

where as is the QCD coupling constant, D^^ is the gluon propagator in the Coulomb gauge 

At: Att ( k^k^\ 

DOO(k) = _^, Z)^.(k) = U^^- - ^ j , D'^ = D^'^ = 0, (6) 

and k = p — q; 7^ and u{p) are the Dirac matrices and spinors 
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with e{p) = y/p'^ + m^. The effective long-range vector vertex is given by 

TM = ^, + ^a,^k\ (8) 

where k is the Pauh interaction constant characterizing the anomalous chromomagnetic 
moment of quarks. Vector and scalar confining potentials in the nonrelativistic limit reduce 
to 

Vv{r) = (1 - e)Ar + B, 

Vsir) = eAr, (9) 

reproducing 

Konf (r) = Vs{r) + Vv{r) = Ar + B, (10) 

where e is the mixing coefficient. 

The expression for the quasipotential for the heavy quarkonia, expanded in f^/c^ without 
retardation corrections to the confining potential, can be found in Ref. [^. The structure 
of the spin-dependent interaction is in agreement with the parameterization of Eichten and 
Feinberg |^T|. All the parameters of our model like quark masses, parameters of the linear 



confining potential A and B, mixing coefficient e and anomalous chromomagnetic quark 
moment k are fixed from the analysis of heavy quarkonium masses (see below Sec. V) and 
radiative decays. The quark masses mi, = 4.88 GeV, rric = 1.55 GeV and the parameters of 
the linear potential A = 0.18 GeV^ and B = —0.16 GeV have usual values of quark models. 
The value of the mixing coefficient of vector and scalar confining potentials e = —1 has 
been determined from the consideration of the heavy quark expansion for the semileptonic 
B D decays and charmonium radiative decays Finally, the universal Pauli 
interaction constant k, = —1 has been fixed from the analysis of the fine splitting of heavy 
quarkonia ^Pj- states ||2^. Note that the long-range magnetic contribution to the potential 
in our model is proportional to (1 + k) and thus vanishes for the chosen value of k = —1. 
In the present paper we will include into consideration the retardation corrections as well 
as one-loop radiative corrections. 



III. GENERAL APPROACH TO ACCOUNTING FOR RETARDATION 
EFFECTS IN THE QQ POTENTIAL 

For the one-gluon exchange part of the qq potential it is quite easy to isolate the retar- 
dation contribution. Indeed due to the vector current conservation (gauge invariance) we 
have the well-known relation on the mass shell 

^Ua{p)Ub{-p)lalbi,Ua{q)Ub{-q) 
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-Ua(p)nf,(-p) 

^0 = ea(p) - ea(q 



« , 1 

k2 A;2 



la ■ lb 



{la ■ k)(7b ■ k) 



k2 



/cq k^; 



p - q 



The left-hand side and the right-hand side of this relation are easily recognized to be in the 
Feynman gauge and the Coulomb gauge, respectively. Now, if the nonrelativistic expansion 
in jvr? is applicable, we can immediately extract the retardation contribution. Namely 
we expand the left-hand side of eq. in 



1 

k2 



A-2 

k4 



and get with needed accuracy 

- na(p)M6(-p) 

In the right-hand side of eq. 
equation 



k2 



1 + 



la-lh 
k2 



Ua(q)u6(-q). 



(12) 



Ur 



p]) one should use the identity following from the Dirac 



ea(q))(eb(q) -eb(p)). 



(p)Mb(-p)(7a ■ k)(7;, ■ k)Ma(q)Mb(-q) 
= 'Ua(p)M6(-p)7°7feMa(q)Mb(-q)(ea(p) 



After defining as a symmetrized product p3| , p4| 



A-2 



(ea(p) -ea(q))(e;,(q) -efe(p)) 



(13) 



and dropping fcg in the denominator we obtain the expression which is identical to eq. 
In this way we obtain the well-known Breit Hamiltonian (the same as in QED |2^) if we 
further expand eq. (|13D in /m? 



/tQ — 



IP - q 



2\2 



AmaTrib 



(14) 



This treatment allows also for the correct Dirac limit in which the retardation contribution 
vanishes when one of the particles becomes infinitely heavy . 



For the confining part of the qq potential the retardation contribution is much more 
indefinite. It is a consequence of our poor knowledge of the confining potential especially in 
what concerns its relativistic properties: the Lorentz structure (scalar, vector, etc.) and the 
dependence on the covariant variables such as k"^ = k^ — k"^. Nevertheless we can perform 
some general considerations and then apply them to a particular case of the linearly rising 
potential. To this end we note that for any nonrelativistic potential l^(— k2) the simplest 
relativistic generalization is to replace it by V{kQ — h"^). 

In the case of the Lorentz-vector confining potential we can use the same approach as 
before even with more general vertices containing the Pauli terms, since the mass-shell vector 
currents are conserved here as well. It is possible to introduce alongside with the "diagonal 
gauge" the so-called "instantaneous gauge" ||2^ which is the generalization of the Coulomb 
gauge. The relation analogous to eq. (|ll|) now looks like (up to the terms of order of /m^) 
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Vv{K - k^)M,(p)M,(-p)r^:r,^u,(q)n,(-q) = Ua{p)u,{-ip)^Vv{-'k')Tyi 

Vy(-k2)r„ • T, + v;(-k2)(r, ■ k)(rfc ■ k)! jM,(q)M,(-q), 



(15) 



where 

Vy{kl - k^) - W(-k^) + klVi.{-\^) 

and as in the case of the one-gluon exchange above we put 



kl = (ea(p) - ea(q))(ef,(q) - eb(p)) 



Amarrib 



(16) 



again with the correct Dirac hmit. 

For the case of the Lorentz-scalar potential we can make the same expansion in /cq, which 
yields 



Vs{K - k') = ^5(-k^) + A:oV^(-k^). 



(17) 



But in this case we have no reasons to fix in the only way (p!3|). The other possibility is 
to take a half sum instead of a symmetrized product, namely to set (see e. g. [p4| , p5| ) 



kn 



:ea(p)-ea(q))' + (eb(q)-eb(p))' 



2 ~'~ ^2 



mf: mt 



(18) 



The Dirac limit is not fulfilled by this choice, but this cannot serve as a decisive argument. 
Thus the most general expression for the energy transfer squared, which incorporates both 
possibilities (0) and (|TB|) has the form 



kl = A(e,(p) - e,(q))(efe(q) - et{p)) + (1 - A)^ [(e„(p) - e^f + (efe(q) - e^ip))' 

where A is the mixing parameter. 

After making expansion in /m? we obtain 



(19) 



(1-A) 



I 1 ^ 



2A 



(k-p)2 + 2(k-p)(k-q) + (k-q)^ 



\m1^ ' ml J rriamh 
Thus as expected k^ ~ 0{p'^/m^) ^ 1. Then the Fourier transform of the potential 

V{kl - k^) ^ V{-\<.^) + klV'{-\if 
with /cq given by eq. (^) can be represented as follows 

2A 



(20) 



rriamb 



X lv(r)p^ + V'(r)-(p-r) 
I r 



w 



(21) 



where {. . denotes the Weyl ordering of operators and 



V{r) 



(27r)^ 



:\/(-k^)e 



2\ „ik-r 



In the case of the one-gluon exchange potential we had A = 1, 



A Anas , . 4: as 



(22) 



(23) 



As for the confining potential we assume it to be a mixture of scalar and vector parts. 
In the nonrelativistic limit we adopt the linearly rising potential 



Voir) = Ar- K)(-k' 



SttA 
(k2p 



(24) 



which we split into scalar and vector parts by introducing the mixing parameter e. The 
possible constant term in Vq has been discussed in [ p!T[| . 



Vo = Vs + Vv] Vs = eVo; Vy = [l - e)Vo. 



(25) 



Hence the retardation contribution (21) from scalar and vector potentials has the form 



:i-A 



2A 



sy 



sy) 



mt mt 



marrib 



{Vsy{r)p' + V!,y{r)^-{p-rr] 



w 



(26) 



where we use the general Ansatz (|T^), (^) for both the scalar and vector potentials for the 
sake of completeness. 

The other spin-independent corrections in our model had been calculated earlier pD|JTT| 



^(1 + 2k) + -\] AVy{r) + {Vy{r)p'} 



w 2 \ rai m 



L){y,(.)p^}„,. (27) 



Adding to the above expression the retardation contributions ( PB] ) and the nonrelativistic 
parts (^) and ( pSj) we obtain the complete spin-independent qq potential: 



Vsi{r) = Vc{r) + Vo{r) + Vyj,{r) + + A[Vc{r) + (1 + 2«:)W] 

8 \ml mi) 



where the velocity-dependent part 

WD(r) = V^^{r) + ^Td W + ^VD(r), 
Vc{r 



2mamb 



1 



4 a, 



1 ) w 



niamb 



x\Vv{r)p^ + V{.{r) 



4 
p ■ r 



2 ~1 



2mamb [ 3 r 
2A 



P 



(P ■ r)' 



w 



ml ml^ 



V 



rriamb 



1-A^ 



w 



mi ml 



(28) 



(29) 
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X < y4r 



Ar 



1 - 



A 



V 



P 



Ay (p ■ r) 



w 



2A 



5 



rriamb 



(1 + A5)p' + (A5-1) 



IP - rj 



eA.c 



Making the natural decomposition 
1 



1 ) w 



WD(r) 



maTJlb 



w 



ml ml^ 



pV,e(r) 



3r 
3r 



+ 



(i-^)(i-^)-^Y 



1 



4 



Ay 
2 

As' 



(1-5)^+5 

^ '2 2 



Ar, 



4 



[(l-£)(l-Av)+£(l-A5)]Ar. 



we obtain from eqs. 

K(r) 

l^de(r) = 7 [(1 - (1 - Ay) - £(1 + \s)] Ar. 
K(r) 

The following simple relations hold: 

He - K = (1 - e)Ar- Vde + K = 
The exact BBP relations |[12| (see also ||2^) in our notations look like 

Kfc - \vtc + \(Vc + Va) = 0, 

2 4 dr 



iP-rj 



(30) 



(31) 



(32) 



(33) 



(in the original version Vhc = — H — and Vde = Vd + |K)- 

The functions ( |3TD identically satisfy the BBP relations ( |33D independently of values of 
the parameters e, Ay, A5 but only with the account of retardation corrections. 

In our model pO| , pA| we have e = —1 and Ay = 1, if we assume further that A5 = 1 |Tl 
then we get 



Hc(r) = + ^Ar; V,{r) 



2a, 1 , 

- -Ar: 

3r 2 ' 
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Vde{r) = -Ar- V^r) = 0. 



(34) 



Our expressions (^Sj) and for purely vector (e = 0) and purely scalar (e = 1) interactions 
and for k = 0, A5 = Ay = 1 coincide with those of Ref. [p5| . 
In the minimal area low (MAL) and flux tube models 



VUr) 



2a, 

3r 
1 



o 

1 



2a, 1 , 

^ - -Ar; 

3r 6 



Kze(r) = --Ar; Ve = --Ar. 
o 



(35) 



To obtain these expressions one should set in relations (|3T|), 

Ay + 2A5 = 1. (36) 
Thus one gets a family of values for Ay and A5. The most natural choice reads as 

(37) 



2 

^=3' 



Ai 



A,. 



0, 



which resembles the Gromes proposal p4[: the symmetrized product for the vector potential 
and the half sum for the scalar potential. But still the Dirac limit is not fulfilled in this case. 
Expression (pSf) for Vsi contains also the term with the Laplacian: 



1/1 1 

+ 



ml ' ml ^ 



A[Vc{r) + {l + 2K)Vv{r) 



In the MAL and some other models these terms look like 28 



o ^ + ^ ^IMr) + Voir) + K(r)] 



m„ 



mt 



and usually it is adopted that 

Lattice simulations [^] suggest that 

AVl'(r) = 



AVa(r) = 0. 



b = 0.8GeV\ 



(38) 



(39) 



(40) 



(41) 



In our model expression (|38D can be recast as follows 

K(r) = (l + 2K)(l-e)V'„(r)-H(r) 

and for the adopted values e = — 1, k = — 1 

AVair) = -3A{Ar) = -6-, 6A = l.lGeV' 



(42) 



(43) 



which is close to the lattice result (^) but differs from the suggestion (|40|) . 
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IV. HEAVY QUARK-ANTIQUARK POTENTIAL WITH THE ACCOUNT OF 
RETARDATION EFFECTS AND ONE LOOP RADIATIVE CORRECTIONS 



At present the static quark-antiquark potential in QCD is known to two loops |T3|,|n[. 
However the velocity dependent and spin- dependent parts are known only to the one-loop 
order [^,0. Thus we limit our analysis to one- loop radiative corrections. The resulting 
heavy quark-antiquark potential can be presented in the form of a sum of spin-independent 
and spin-dependent parts. For the spin- independent part using the relations (pSD, (^) with 
Ay = 1 and including one-loop radiative corrections in MS renormalization scheme p!5|JT6 
we get 



V^si(r) 



+ 



r 



3 
1 



^ 4 /3oag(/i^)ln(/ir) 
3 2tx 



A 



mt 



+ 



2mamh \ [ 

3 2n 
1-e 
2mamb 



P 



4 ay 
3 r 



4 ay (/i^) 4 Poa'j, (/i^) ln(/ir ) 



27r 



+ {l-e){l + 2K)Ar 



r 

e 
4 



ml 



+ 
1 

2 

mi 



J J w 

(p-r)V ln(/xr) l\] 



Ar 



P 



IP-rj 



w 



+ 



mt 



m 



4 \ml 



1 1 



b , 
+ 



rrianib 



Ar 



5p^ 



where 



1 + 



4 2 



TT 



(44) 



(45) 
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ai = rif, 

3 9 



II- -n. 



Here is a number of flavours and /i is a renormalization scale. 

For the dependence of the QCD coupling constant as(/i^) on the renormalization point 
/i^ we use the leading order result 



All 



/3oln(^VA2) 

Comparing this expression for Vsi with the decomposition (0) we find 
2ay(/i^) 2 /?oas(/i^) ln(/ir) /l— ^ ^^s^ 



(46) 



Hc(r) = - 



3 r 3 27r 

2 «vV) _ 2^o«|(/i^ 

3 r 3 27r 



r 

ln(/ir) 



2 

'1 



I + 5, 
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Ve{r) 



'l + Xs)Ar + B, 



--(l-A5)Ar. 



(47) 



It is easy to check that the BBP relations are exactly satisfied. 

The spin-dependent part of the quark-antiquark potential for equal quark masses (m^ = 
m}y = m) with the inclusion of radiative corrections pr5| , p!7[| can be presented in our model 
as follows: 



SD 



a L- S + 6 

f4as(/i 



— (Sa ■ r)(Sb ■ r) 



2m2 j' 



1 + 



TT 



1 



— Uf h 7k 

18 ^ 36 ' 



+ (^^ - 2^ In(mr) ^ - ^ + A{1 + k){1 - e)^ 



1 

J ^4 



3 



1 + 



In (mr) 



TC 



25 
12 



+ a + K) 



+ lE 



2 



r 



87ras(/i^ 



3m^ [ 



1 + 



TT 



23 5 3 

Uf m 2 

12 18 ^ 4 



vr 



+ (l + /«) 



87r \ 

r 



'ln(/i/m) 



1 

^ V12 



1 1 

Uf - 



16 



+ ^lnii 
2 m 



^Inii 
2 m 



I In(mr) + 7^;' 



(48) 



(49) 



(50) 



(51) 



where L is the orbital momentum and ^a,h-, S = Sq + Sf, are the spin momenta. 

The correct description of the fine structure of the heavy quarkonium mass spectrum 
requires the vanishing of the vector confinement contribution. This can be achieved by 
setting 1 + K = 0, i.e. the total long-range quark chromomagnetic moment equals to zero, 
which is in accord with the fiux tube and minimal area ||30| , |28[| models. One can see from 
Eq. (^) that for the spin-dependent part of the potential this conjecture is equivalent to 
the assumption about the scalar structure of confinement interaction 



V. HEAVY QUARKONIUM MASS SPECTRA 

Now we can calculate the mass spectra of heavy quarkonia with the account of all rel- 
ativistic corrections (including retardation effects) of order f^/c^ and one-loop radiative 
corrections. For this purpose we substitute the quasipotential which is a sum of the spin- 
independent (^iD and spin-dependent (|48|) parts into the quasipotential equation (|l]). Then 
we multiply the resulting expression from the left by the quasipotential wave function of a 
bound state and integrate with respect to the relative momentum. Taking into account the 
accuracy of the calculations, we can use for the resulting matrix elements the wave functions 
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of Eq. (|l|) with the static potential [| 



VMr) = --^-^ + Ar + B. 
6 r 

As a result we obtain the mass formula {nia = rrib = m) 



b\M) 



W+{a){L-S) + {b){ 



;{Sa ■ r)(Sfe ■ r) - {Sa ■ Sfe 



) + (c)(S,-S,), 



(52) 



(53) 



where 



W = (Vsi) + 
1 



( 



-(S,T)(Sfe-r)-(S„-S 



L • S) = -(J(J + 1) - L{L + 1) - ^(^ + 1)), 

6((L • S))2 + 3(L ■ S) - 25(5 + 1)L{L + 1] 



) 



2(2L- l)(2L + 3) 
(S, -86) = ^ (5(5 + 1) - ^) , S = + S,, 



and (a), (6), (c) are the appropriate averages over radial wave functions of Eqs. (^Dp-(^Tp. 
We use the usual notations for heavy quarkonia classification: where n is a radial 

quantum number, L is the angular momentum, 5 = 0, 1 is the total spin, and J = L — 
S, L, L + S is the total angular momentum (J = L + S). 

The first term on the right-hand side of the mass formula (^3|) contains all spin- 
independent contributions, the second term describes the spin-orbit interaction, the third 
term is responsible for the tensor interaction, while the last term gives the spin-spin inter- 
action. 

To proceed further we need to discuss the parameters of our model. There is the following 
set of parameters: the quark masses (m^ and rric), the QCD constant A and renormalization 
point /i (see Eqs. (^61), (^41), (^S])) in the short-range part of the QQ potential, the slope A 
and intercept B of the linear confining potential (|T3), the mixing coefficient e (y), the long- 
range anomalous chromomagnetic moment k of the quark (|^), and the mixing parameter 
Xs in the retardation correction for the scalar confining potential (^). As it was already 
discussed in Sec. II, we can fix the values of the parameters e = —1 and n = —1 from the 
consideration of radiative decays and comparison of the heavy quark expansion in our 



model |2^j3^ with the predictions of the heavy quark effective theory. We fix the slope of 
the linear confining potential A = 0.18 GeV^ which is a rather adopted value. In order to 
reduce the number of independent parameters we assume that the renormalization scale /i in 
the strong coupling constant Q;s(/i^) is equal to the quark mass. | We also varied the quark 



^This static potential includes also some radiative corrections |T^. The remaining radiative cor- 
rection term with logarithm in (]4^), also not vanishing in the static limit, is treated perturbatively. 

^Our numerical analysis showed that this is a good approximation, since the variation of fi does 
not increase considerably the quality of the mass spectrum fit. 
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masses in a reasonable range for the constituent quark masses. The numerical analysis and 
comparison with experimental data lead to the following values of our model parameters: 



= 1.55 GeV, = 4.88 GeV, A = 0.178 GeV, 
A = 0.18 GeV^ B = -0.16 GeV, fi = mg {Q = c,b), 
e = —1, K = —1, \q = 0. 



The quark masses mc^b have usual values for constituent quark models and coincide with 
those chosen in our previous analysis ||2^ (see Sec. II). The above value of the retardation 
parameter for the scalar confining potential coincides with the minimal area low and 



flux tube models |28|, with lattice results |^ and Gromes suggestion [gj]. The found 
value for the QCD parameter A gives the following values for the strong coupling constants 
asiml) ^ 0.32 and as{ml) ^ 0.22. 

The results of our numerical calculations of the mass spectra of charmonium and bot- 
tomonium are presented in Tables | and ||. We see that the calculated masses agree with 
experimental values within few MeV and this difference is compatible with the estimates 
of the higher order corrections in f^/c^ and a^- The model reproduces correctly both the 
positions of the centres of gravity of the levels and their fine and hyperfine splitting. Note 
that the good agreement of the calculated mass spectra with experimental data is achieved 
by systematic accounting for all relativistic corrections (including retardation corrections) 
of order v'^/c^, both spin-dependent and spin- independent ones, while in most of potential 
models only the spin-dependent corrections are included. 

The calculated mass spectra of charmonium and bottomonium are close to the results 
of our previous calculation ||2^ where retardation effects in the confining potential and 
radiative corrections to the one-gluon exchange potential were not taken into account. Both 
calculations give close values for the experimentally measured states as well as for the yet 
unobserved ones. The inclusion of radiative corrections allowed to get better results for the 
fine splittings of quarkonium states. Thus we can conclude from this comparison that the 
inclusion of retardation effects and spin-independent one-loop radiative corrections resulted 
only in the slight shift (~ 10%)in the value of the QCD parameter A and an approximately 
two- fold decrease of the constant B. | Such changes of parameters almost do not influence 
the wave functions. As a result the decay matrix elements involving heavy quarkonium 
states remain mostly unchanged. We plot the reduced radial wave functions u{r) = rR{r) 
for charmonium and bottomonium in Figs. |I| and ^ 



Note that in Ref. [20| we included this constant both in vector and scalar parts, while the present 
analysis indicates that the better fit can be obtained if the constant B is included only in the vector 
part (|). 

^The changes in decay matrix elements are of the same order of magnitude as the contributions 
of the higher order relativistic and radiative corrections. 
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VI. CONCLUSIONS 



In this paper we have considered the heavy quarkonium spectroscopy in the framework 
of the relativistic quark modeL Both relativistic corrections of order v"^ /(? and one-loop 
radiative corrections to the short-range potential have been included into the calculation. 
Special attention has been devoted to the role and the structure of retardation corrections 
to the confining interaction. Our general analysis of the retardation effects has shown 
that we have a good theoretical motivation to fix the form of retardation contributions to 
the vector potential in the form ([ISD which corresponds to the parameter Ay = 1 in the 
generalized expression (P^l). On the contrary, the structure of the retardation contribution 
to the scalar potential is less restricted from general analysis. This means that it is not 
possible to fix the value of A5 in ( PBD on general grounds. Our numerical analysis has 
shown that the value of As' = is preferable. Thus for the energy transfer squared we 
have the symmetrized product (^) for the vector potential and a half sum (|T^) for the 
scalar potential, in agreement with lattice calculations |29| and minimal area law and flux 



tube models [^. The found structure of the spin- independent interaction (^4]) with the 
account of retardation contributions satisfies the BBP |12[ relations (0), which follow from 
the Lorentz invar iance of the Wilson loop. 

In our calculations we have used the heavy quark-antiquark interaction potential with 
the complete account of all relativistic corrections of order f^/c^ and one-loop radiative 
corrections both for the spin-independent and spin-dependent parts. The inclusion of these 
corrections allowed to fit correctly the position of the centres of gravity of the heavy quarko- 
nium levels as well as their fine and hyperfine splittings. Moreover, the account for radiative 
corrections results in a better description of level splittings. The values of the main param- 
eters of our quark model such as the slope of the confining linear potential A = 0.18 GeV^, 
the mixing coefficient e = —1 oi scalar and vector confining potentials and the long-range 
anomalous chromo magnetic quark moment k = — 1 used in the present analysis are kept the 
same as they were fixed from the previous consideration of radiative decays |H] and the heavy 



quark expansion ||22|,^. The value oi e = —1 implies that the confining quark-antiquark 
potential in heavy mesons has predominantly a Lorentz-vector structure, while the scalar 
potential is anticonfining and helps to reproduce the initial nonrelativistic potential. On the 
other hand, the value of k = — 1 supports the conjecture that the long-range confining forces 
are dominated by chromoelectric interaction and that the chromomagnetic interaction van- 



ishes, which is in accord with the dual superconductivity picture |^ and flux tube model 



The presented results for the charmonium and bottomonium mass spectra agree well 
with the available experimental data. It is of great interest to consider the predictions for 
the masses of the ^Sq and D levels of bottomonium, which have not yet been observed 
experimentally. The difficulty of their experimental observation is that these states (except 
^Di) cannot be produced in e^e~ collisions, since their quantum numbers are not the same 
as the quantum numbers of the photon. Therefore, in search for these states one must 
investigate decay processes of vector i^Si) levels. We discussed the possibility of observation 
of these states in radiative decays in Ref. |^. Note that the small value predicted for the 
hyperfine splitting M(T) — M{rii)) = 60 MeV leads to difficulties in observation of the rjb 
state. 
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Recently it was argued that the account of relativistic kinematics substantially modi- 
fies the description of the charmonium fine structure and, in particular, leads to considerably 
larger values of the 2^Pj splittings than in the nonrelativistic limit. Both our previous calcu- 
lation pOf and the present one confirm this observation. Our prediction for the charmonium 
2^Po mass lies close to the prediction of Ref. |^ and slightly lower than the DD* threshold. 
However, the fact that this state is above DD and close to DD* thresholds makes threshold 
effects very important and can considerably infiuence the quark model prediction. 
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TABLES 



TABLE I. Charmonium mass spectrum. 



State (n(2S+i)Lj) 



Particle 



Theory 



Experiment [1| 



Experiment |31] 



Vc 



2.979 
3.096 



2.9798 
3.09688 



2.9758 



XcO 
Xcl 
Xc2 



3.424 
3.510 
3.556 



3.4173 
3.51053 
3.55617 



3.4141 



2^50 
23^1 



V'c 



3.583 
3.686 



3.594 
3.686 



3.798 
3.813 
3.815 



3.7699* 



23Po 
2^ Pi 

2^P2 



X'cO 
X'cl 
Xc2 



3.854 
3.929 
3.972 



3^50 
3351 

2^D2 
2^D. 



Vc 



3.991 
4.C 



4.194 
4.215 
4.223 



4.040* 
4.159* 



* Mixture of S and D states 



18 



TABLE II. Bottomonium mass spectrum. 



State (n(2S+i)Lj) 



Particle 



Theory 



Experiment [I| 



Experiment |32] 



T 



9.400 
9.460 



9.46037 



Xbo 
Xbi 

Xb2 



9.864 
9.892 
9.912 



9.8598 
9.8919 
9.9132 



9.8630 
9.8945 
9.9125 



2^ So 
2^Si 



9.990 
10.020 



10.023 



10.151 
10.157 
10.160 



23Po 
2^ Pi 

2^P2 



x'bo 
X'bi 

X'b2 



10.232 
10.253 
10.267 



10.232 
10.2552 
10.2685 



3^50 
3351 



10.328 
10.355 



10.3553 



2-Wi 
2^D2 
23L», 



10.441 
10.446 
10.450 



3^Po 
33 Pi 

33P2 



Xbo 
X'li 

X'l2 



10.498 
10.516 
10.529 



4^50 
4351 



10.578 
10.604 



10.580 
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FIGURES 




r-GeV 

FIG. 2. The same as in Fig. 1 for bottomonium and long-short-dashed hne for 3S state. 
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